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We present a simplified computational scheme in order to calculate the effects of electron correla- 
tions on the energy bands of diamond and silicon. By adopting a quasiparticle picture we compute 
first the relaxation and polarization effects around an electron set into a conduction band Wannier 
orbital. This is done by allowing the valence orbitals to relax within a self-consistent field (SCF) 
calculation. The diagonal matrix element of the Hamiltonian leads to a shift of the center of gravity 
of the conduction band while the off-diagonal matrix elements result in a small reduction of the 
conduction-electron band width. This calculation is supplemented by the computation of the loss 
of ground state correlations due to the blocked Wannier orbital into which the added electron has 
been placed. The same procedure applies to the removal of an electron, i.e., to the valence bands. 
But the latter have been calculated previously in some detail and previous results are used in order 
to estimate the energy gap in the two materials. The numerical data reported here shows that the 
methods works, in principle, but that also some extension of the scheme is necessary to obtain fully 
satisfactory results. 



Dedicated to J. -P. Malrieu on the occasion of his 60th birthday. 



I. INTRODUCTION 



Reliable energy-band calculations of solids remain a 
major topic of modern condensed-matter theory. The 
problem is a very old one and has been treated on vari- 
ous levels of accuracy and sophistication. Originally the 
calculations were based on a one-electron theory, i.e., 
many-body effects were simply neglected or indirectly ac- 
counted for 1 . A good choice of model potentials was cru- 
cial for the determination of energy bands, e.g., of semi- 
conductors or metals. After the development of density 
functional theory (DFT) the field obtained an immense 
impact. Although the theory is designed for ground-state 
properties it has been also used for the determination 
of energy bands. Within the local density approxima- 
tion (LDA) density functional theory could be cast into 
the form of the Kohn-Sham equation 2 containing a self- 
consistent potential which is free of model-dependent as- 
sumptions. The eigenvalues of that equation are sim- 
ply interpreted as representing the energy bands of the 
system. The computational efforts were significantly re- 
duced and simplified by the development of linearized 
methods 3 . The results agreed often surprisingly well 
with experiments indicating partial error cancellations. 
In cases where the agreement is poor it is difficult to 
make improvements which are free of arbitrariness. For 
example, for strongly correlated electrons the LDA+C7 
method has been applied 4 where the intra-atomic corre- 
lations are described by a Hubbard interaction term with 
an on-site Coulomb integral U. Thereby the problem of 



double counting of the interactions must be carefully con- 
sidered. More recently the LDA+U has been combined 
with the dynamical mean-field theory (DMFT) 5 , an ex- 
tension of the coherent potential approximation (CPA) to 
include the effects of dynamical fluctuations (dynamical 
CPA) 6 , or the so-called GW method 7 - 8 . 

One of the shortcomings of the LDA when applied to 
semiconductors or insulators is the underestimation of 
the size of the gap between the valence and conduction 
bands. It is of no surprise that the LDA fails here because 
the correlation cloud around an electron set into the con- 
duction band differs strongly from that of electrons in the 
charge-neutral ground state. While the former includes a 
long-ranged polarization cloud the latter consists of van 
der Waals fluctuations. But those two different physical 
situations cannot be described by a single self-consistent 
potential depending on the density only as assumed in 
the LDA. Therefore one would like to apply quantum- 
chemical methods in order to calculate energy bands of 
semiconductors. They allow for controlled approxima- 
tions and can be done at different levels of accuracy. 

Calculations based on quantum-chemical methods re- 
quire a self-consistent field (SCF) or Hartree-Fock (HF) 
energy-band calculation as a starting point. With the 
CRYSTAL program package 9 such calculations can be 
done. Furthermore, because of the local character of 
the correlation hole around an electron one also needs 
to know the Wannier functions for a given semiconduc- 
tor. Correlation effects arise from exciting electrons out 
of those localized states. Wannier functions can be ob- 
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taincd either from the CRYSTAL program package 9 or 
from an alternative code called WANNIER 10 . While it 
is straight forward to determine Wannier orbitals for va- 
lence bands it is much more difficult to do the same for 
conduction bands, due to their weaker localization and 
the frequently encountered mixing of the low-lying con- 
duction bands with energetically higher-lying bands. But 
with a proper band disentanglement procedure this can 
be achieved too 11 . It would be preferable to excite elec- 
trons out of local non-orthogonal orbitals 12 implying the 
use of a different set of orbitals for a SCF and a cor- 
relation calculation. But in that case one cannot use 
standard quantum chemical programs but rather has to 
develop new ones. 

In order to calculate energy bands for semiconduc- 
tors with quantum chemical methods one can proceed 
along two different routes. One starts from the quasi- 
particle picture 13,14 . The real part of the excitation en- 
ergy is assumed to be much larger than the imaginary 
part so that the life time of the excitation may be ne- 
glected. In that case one identifies the most important 
microscopic processes which generate the correlation hole 
around an electron or hole added to the system and di- 
agonalizes the Hamiltonian in that particular subspace. 
The latter includes the additional electron or hole (pro- 
jection method). A related method is based on an effec- 
tive Hamiltonian and has been exploited by Malrieu and 
coworkers 15,16 and others 17 . The second route consists 
of calculating Green's functions 18-22 , the poles of which 
determine the excitation spectrum. In both cases an in- 
cremental scheme 23 can be applied when the correlation 
hole is constructed around the added particle or hole. A 
number of high-quality calculations have been performed 
for the valence energy bands, e.g., of diamond and sil- 
icon 14 ' 24 ' 25 . However, one is still not able to calculate 
the conduction bands (and thus the energy gap) of these 
materials in the same way, even with our new embedding 
tool 26,27 Convergence of the approach to an accuracy bet- 
ter than 0.5 eV seems presently hard to achieve. 

Here we want to suggest an alternative scheme for en- 
ergy bands which is based on the quasiparticle approach. 
It uses a division of the correlation effects caused by the 
added particle (i.e., electron or hole) into two parts. One 
concerns polarization and relaxation around the particle 
and another the loss of ground-state correlations associ- 
ated with its presence. A splitting of this form was stud- 
ied previously for a simplified model which allowed for a 
nearly analytic treatment of those correlation contribu- 
tions to the energy bands 28 . Here we want to demon- 
strate that it can be used also for more accurate results 
based on an ab initio Hamiltonian for a respectable ba- 
sis set. Calculations of this kind seem to be much easier 
than previously used computational schemes for dealing 
with the Hamiltonian. 

The numerical results which we shall present here re- 
quire further work and improvements. That is related to 
the rather diffuse form of the Wannier functions for con- 
duction electrons in the covalent systems diamond and 



silicon. But it is also demonstrated that the improve- 
ments are straight forward. 

II. DESCRIPTION OF THE METHOD 

We assume that the Hamiltonian is written for a ba- 
sis set consisting of Gauss- type orbitals (GTO's). Fur- 
thermore, we assume that the Hartree-Fock (HF) or self- 
consistent field ground-state wave function |$scf) of the 
TV-electron system has been determined. For the el- 
emental semiconductors we are interested in, this can 
be achieved by using the program package CRYSTAL 9 . 
When an electron with momentum k and spin a is added 
to the energy band v the (A+I)-particle state is 

|<e +1 (k))=c+ J* SCF ) • (I) 

Similarly, when an electron is removed from the system 
the state is 

|^- 1 (k))=c k , ff |$ SCF ) • (2) 

Next wc want to express the Bloch states as superposi- 
tion of Wannier states. From the occupied HF orbitals 
one can construct Wannier orbitals without problems, 
e.g., by using the Wannier-Boys localization method 29 
implemented in CRYSTAL 9 . For the virtual space this 
is more difficult. The required procedure is described in 
Rcf. 11 to which we refer. In the following we assume 
that those orbitals have been determined and that the 
corresponding creation (destruction) operator are c^ n(T 
(cRncr)- Here R denotes the lattice vector of the unit cell 
at which the Wannier orbital is centered while n is an 
intracell index. One can imagine the Wannier orbitals as 
representing the bonding and antibonding orbitals in the 
diamond structure of the elemental semiconductors. Wc 
write 

| < ff +1 (k)) = -j= ]rcw(kK kR c+„ a | $scf) (3) 

' c Rn 

and similarly for |<J>^ T ~ 1 (k)). Here N c is the number of 
unit cells. The conduction-band energies in Hartree-Fock 
approximation can be expressed in terms of local matrix 
elements as 

47 = EE«-( k )^( k ) eJkR «' ( 4 ) 

R nn' 

with matrix elements 

#R-R,rm' = (*SCF | CR„' CTJ ffc+ mJ | $ S Cf) 

~<5rR 5 nn > Eq CF . (5) 

Here Eq CF is the Hartree-Fock ground-state energy. 
When the electronic correlations are included, the 
Hartree-Fock excitations become quasiparticle excita- 
tions. For relatively weakly correlated electron systems 
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one may assume that there is still a one-to-one correspon- 
dence of the Hartree-Fock and quasiparticle excitations. 
In this case the energy bands describe the dispersion of 
the quasiparticlcs. Life-time effects are neglected. Thus 
we write for the correlated (7V+l)-electron state 



<+ 1 (k))=n|<+ 1 (k)) 



5> ra (k)e ikR fic+ n( , 

Rti 



$scf) (6) 



where the wave- or Moller operator Q transforms 
|$^ +1 ( k )> into l*^r +1 ( k )>- Thc correlated conduction 
bands are given by 



e (h) _ (^ +1 (k)|»+g»|^+i(k)) 

vaK} (efkjisi+fiiefk)) 



(7) 



where Eq is the exact, i.e., correlated ground-state en- 
ergy. In terms of cumulants this can be rewritten as 

e„ CT (k) = (H | Q) ki , CT - E (8) 

where the notation 

(A | B) kva = ($^ +1 (k) | A+B | <C +1 (k)) c (9) 

has been used with the index c indicating that the cumu- 
lant has to be taken. For cumulants and their properties 
we refer to Ref. 30 where also the original literature is 
cited. The operator fi is the cumulant wave operator, i.e., 
it is defined only in connection with the form (9). It is 
customary to write = 1 + 5 where S is a scattering ma- 
trix. By decomposing the Hamiltonian H = Hscf + ^rcs 
we obtain for the correlation correction of the conduction 
bands 

e w (k) =e^(k) + Ac w (k) 
Ae w (k) = ^^« ra (kK n ,(k)e' kR 

R nn' 

(c+ v | h ics s CR „ ff ) . (10) 

The round brackets refer to the cumulants 



(A | B) = ($scf | A+B | $ SCF ) C 



(11) 



Note that the cumulant (c+ n , a \ H les S cj^ ncr ) contains a 
term -E CO rr = Eq — Eq CF in form of a subtraction which 
describes the correlation energy of thc ground state of thc 
TV-electron system. It changes when differently localized 
Wannier functions are used. Yet, since Wannier func- 
tions provide a complete basis, the shifts Ae va remain 
the same after summation of all relevant contributions 
according to (10). 

For systems like the elemental semiconductors it suf- 
fices to include in S one- and two-particle excitations. 
Thereby, thc annihilation operators CR n(T refer to the oc- 
cupied Wannier orbitals of thc (A+l)-clcctron system 
while the number of the local creation operators a^ T(T 



depends on the size of the basis set. The latter are, e.g., 
the ones for GTO's with the GTO being centered at site 
(or bond) r of the unit cell R. Note that due to the use of 
cumulants only connected operator products contribute. 
Thus the following general ansatz for the scattering ma- 
trix S holds 



ma 

R1R2 Tn a a' 

1 R.2 R-3 R-4 

R1R2R3R4TT'™' uc' 
a Ri r a a R 2 r ' a ' C R-3 ncr ' C R 4 n ' a 



(12) 



It is advantageous to decompose 5* into a part which com- 
mutes with the added electron c^ ncr and a part which 
does not, i.e., 



S = + S v (13) 



with 



[S„cL] _ ¥= and [S v , c+ na ] _=0 . (14) 

The operator SV describes the relaxation and polariza- 
tion around thc electron generated by ci while S v de- 
scribes the loss of ground-state correlations due to pres- 
ence of the extra electron. As pointed out before, thc 
latter results from blocking thc orbital occupied by the 
added electron for excitations out of the ground state. 
Thus 

( C 0n>* I #rcs S C+ na ) = ($ S CF | C 0n , a H Ies S w C+ n „ | $ SC f) C 
+ ($SCF I Con^H^C^Sr, | $ S Cf) C - 

(15) 

The first term, i.e., the polarization and relaxation energy 
is approximated by a SCF calculation for the (N+l)- 
electron system with an electron frozen in Wannier or- 
bital (Una). This can be done by making use of the 
program package MOLPRO . The energy contribution 
of the second term is computed by determining those 
excitations in a calculation of the ground-state correla- 
tion energy which are blocked by the creation c^ ncr of 
an electron. Their contribution to thc energy can be 
extracted from a ground-state calculation. The details 
of this procedure are described in thc following subsec- 
tion. Needless to say that the same prescription applies 
for calculating thc valence band, i.e., when an electron is 
removed from the iV-particle ground state. 

Although the scheme proposed here consists of sev- 
eral individual steps and may seem cumbersome from a 
technical point of view, the present approach is concep- 
tionally much simpler than other ab initio methods for 
band structures 14 ' 21 ' 22 ' 24 ' 25 ' 32 . It is designed such that 
it only requires Hartree-Fock type calculations and some 
small configuration interaction (CI) correlation calcula- 
tions to estimate the correlation effects on the energy 
bands. Both type of calculations can be performed rou- 
tinely and without much computational effort by stan- 
dard quantum chemical program packages. 
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III. COMPUTATIONS AND RESULTS 

In a series of papers 14 ' 24 ' 25 , we presented results for the 
influence of correlation effects on the valence band struc- 
ture of diamond and silicon. These results were obtained 
using a local Hamiltonian formalism in finite cluster cal- 
culations for the neutral iV-particle state and cationic 
(TV-l)-particle hole states. Unfortunately, the structure 
of the conduction band is not accessible this way, since an 
additional electron in negatively charged (7V+l)-particlc 
states of finite C and Si clusters is not bound. Therefore, 
only rough estimates of the band gaps in diamond and 
silicon could be given in Ref. 25, assuming that the ener- 
getic effect of the polarization cloud around a conduction 
band electron would be approximately the same as that 
around a valence-band hole. 

In the meantime, we adapted our local Hamiltonian 
formalism to the treatment of conduction bands, by ex- 
tracting the localized conduction band Wannier functions 
(WF) from periodic Hartree-Fock (HF) calculations 29 ' 33 . 
The same calculations are used to generate embedding 
potentials which describe the influence of the infinite 
bulk surroundings on finite cluster models 26 . Results for 
the band structure of diamond 34 and hydrogen fluoride 
chains 35 with this approach are under preparation in our 
laboratory. 

As an alternative, we present here results for the di- 
amond and silicon conduction bands using the simpli- 
fied method described in the previous section. Thereby, 
instead of generating local matrix elements of the sys- 
tem's Hamiltonian between localized multi-determinant 
conduction band states, we start from the simple physical 
picture of a polarization cloud forming around an added 
conduction-band electron and calculate the (static) po- 
larization and relaxation effects around this electron 
when it is placed in a frozen conduction band WF. This 
yields localized quasi-particle states which we can use as 
a basis for the final diagonalization of the Hamiltonian. 
This in turn should lead to a realistic quantitative de- 
scription of correlation effects on the conduction bands. 
The advantage over the local Hamiltonian approach lies 
in the facts that by freezing the conduction band WF, 
we do not encounter stability problems, and that for the 
determination of the static (rather than dynamic) po- 
larization cloud self-consistent field (SCF) calculations 
(rather than correlated calculations) are sufficient. It is 
clear, on the other hand, that the computational sim- 
plification cannot be achieved without introducing ad- 
ditional approximations: the difference between static 
and dynamic relaxation/polarization effects is neglected 
and, more importantly, correlation effects beyond relax- 
ation/polarization are also neglected at this stage. The 
most important of the latter effects is the loss of ground- 
state correlation when introducing an electron in the con- 
duction band, i.e., the exclusion effect due to the non- 
availability of excitations from occupied orbitals of the 
neutral (TV-electron) system into the conduction band 



Wannier spin orbital occupied in the (iV+l)-electron sys- 
tem. However, the leading term of this exclusion effect 
on the diagonal matrix clement can easily be recovered 
in a small configuration-interaction (CI) calculation, re- 
stricted to excitations from the occupied SCF orbitals of 
the closed-shell A-electron system into a single localized 
conduction band WF, without abandoning the advan- 
tages of computational simplicity. 

Following the above strategy, we performed the follow- 
ing detailed computational steps: 

1) We started from periodic HF calculations for 
bulk diamond and silicon. The CRYSTAL200x 
program 33 was used for that purpose, a pre-cursor 
of the most recent public version of the CRYSTAL 
program package 9 which already allowed for local- 
ization. The basis sets used were of valence-triplc- 
zeta quality, with a single d polarization function, 
optimized for diamond 36 ' 37 and silicon 38 . In the 
case of Si, the [Is2s2p] core was simulated by an 
effective core potential 38 . These calculations pro- 
vide the HF band structure of the conduction (and 
valence) bands. 

2) The Wannicr-Boys algorithm 29,39 is used to opti- 
mize the unitary band mixing matrix of the multi- 
band Wannier transformation. It is applied sep- 
arately to the valence band and conduction band 
Bloch functions to generate localized occupied and 
virtual Wannier functions. In the case of Si, the 
conduction band states had to be disentangled from 
those of the higher-lying unoccupied bands be- 
fore they could be localized 11 . The mean width 
\J ((r — (r)) 2 ) of the resulting virtual Wannier func- 
tions is found to be in the order of 1.5-2 times the 
bond length in the elemental semiconductor. 

3) The data of steps 1 and 2 were transferred via our 
CRYSTAL-MOLPRO interface 27 to the molecular 
program MOLPRO 31,40 and are used as input for 
generating embedded C 8 ( 18 ) and Si 8 ( 18 ) clusters 34 . 
Here, the central C-C or Si-Si bond and its six 
nearest-neighbor bonds inside the C 8 (Si 8 ) kernel 
are treated as active, while the basis functions from 
the surrounding shell of 18 C (Si) atoms are only 
added to assist in the representation of the local- 
ized occupied and virtual orbitals of the kernel. 
These C26 and Si26 support clusters are sufficient to 
provide good representations of the localized Wan- 
nier functions inside the kernel. Upon projection 
onto the support clusters the norms of the virtual 
Wannier functions reduce by 1-2% only. For the 
generation of the local virtual basis functions a 
compact variant 26 of the projected atomic orbital 
construction 41 is used. The Hartree-Fock potential 
of all non-active occupied C-C/Si-Si orbitals on the 
cluster and the surrounding solid (together with the 
nuclei there) is incorporated in the embedding po- 
tential. 
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4) The HF reference data for matrix elements between 
localized (7V+l)-particle states can now be gener- 
ated according to 



(*»+ 1 \H\*? + i ) - <$SCF|ff|$SCF)4 



<N+1 



(16) 



where $scf is the HF wave function of the neutral 
^8(18) or Si8(i8) cluster, and <J>f +1 is the same wave 
function with a localized (anti-bonding) conduction 
band WF (from step 2) added at bond i. The lat- 
ter describes the local HF configurations cj^ no .<£>sCF 
introduced in Eq. (3). 

5) To estimate the correlation effects for the diago- 
nal matrix element (i = j), we first replace $^ +1 
in Eq. (16) by $| JV+1 ^ where the valence orbitals 
are allowed to relax and polarize due to the pres- 
ence of the (fixed) conduction band WF of the 
(A r +l)-particle system. This single-determinant 
wave function serves as a first approximation to the 
correlated counter parts ^c^ II(T $scf from Eq. (6). 

6) To include the leading term to the loss of ground 
state correlation for the diagonal matrix element, 
we further replace l>f +1 by and $ S cf by 
V N in Eq. (16), where ^f +1 and * w are full CI 
wave functions in the space of the active valence 
band orbitals (from step 4 for the N-particle sys- 
tem, and step 5 for the (N+l)-particle system) plus 
the localized conduction band WF at site i. (We 
also tried including additional excitation channels 
with one electron in one of the conduction-band 
WF j adjacent to i, but the effect on the numerical 
results to be described below was negligible [<0.05 
eV].) 

7) In order to get the leading term for the correla- 
tion effects on the off-diagonal matrix element be- 
tween nearest-neighbor bond sites i ^ j, we focus 
on the relaxation and polarization effect and ap- 
proximately evaluate 



-;jv+i 



JV+li 



S AT+1\ 



The resulting SCF and correlation contributions are 
listed in Table 1. The employed Wannier functions are 
rather compact, and the matrix elements of the local 
Hamiltonian should decay rapidly with increasing dis- 
tance of the bond pairs, such that the diagonal and 
nearest-neighbor off-diagonal terms should suffice to re- 
cover most of the correlation effects. Switching to differ- 
ent localized Wannier functions would alter these contri- 
butions slightly, but after diagonalization of the result- 
ing Hamiltonian matrix these modifications cancel to a 
large extent. At present, we do not have any alterna- 
tive localization scheme for Wannier functions at hand 
to corroborate this numerically. 

The above diagonal correlation contributions lead to a 
shift of the center of gravity of the band and sum up to 
-1.08 eV for diamond and -0.85 eV for silicon. These en- 
ergies are significantly smaller in magnitude than the cor- 
responding values found in Rcf. 25 for the valence band 
(relaxation/polarization: 4.6 eV for C, 3.4 eV for Si; loss 
of ground-state correlation: -0.7 eV for C, -0.6 eV for 
Si). This finding is in strong contrast to an assumption 
put forward in Ref. 25, namely that the correlation ef- 
fects caused by conduction electrons and valence holes are 
symmetric and presumably has to be attributed to the 
fact, that the conduction band WFs are much less local- 
ized than the valence band WFs. Still, the assumption 
is bound to hold for long-range polarization effects, since 
the polarizing charges are of the same magnitude (±e), 
and if we extract the polarization contributions beyond 
nearest-neighbor bonds (not calculated in the present pa- 
per) from the valence-band calculations of Ref. 25 (3.5 eV 
for C, and 2.4 eV for Si, without two- and three-body con- 
tributions), we can estimate the total correlation- induced 
shift of the conduction band center of gravity to -4.6 eV 
for C, and -3.3 eV for Si. A much simpler estimate of 
the long-range polarization contribution to the shift of 
the energy bands can be obtained from the classical po- 
larization energy of a dielectric medium outside a sphere 
of radius i? cu t (see Ref. 30) 



+2($ 



N+V 



H\$ 



N+l 



- $ 



N+l 



) (17) 



eo-1 e 2 
2eo i?cut 



(19) 



with eo being the dielectric constant of the medium. Us- 
ing the experimental values for the dielectric constant 



exploiting the fact that $j and ^- are equiv- 
alent in diamond lattices. To this end the polarized 
active valence orbitals (p n entering the determinant 
&f +1 are expanded into the unpolarized ones (ip n ) 

defining &f +1 plus orthogonal relaxation contribu- 
tions u>k in the virtual space 

<fin = V 1 - C l <Pn + C n 0J n , (18) 

and is replaced by a CI wave function with 

single excitations from the <p n into the u) n . 



TABLE I. Contributions to the diagonal and off-diagonal 
Hamiltonian matrix elements (in eV) between localized con- 
duction band states in bulk diamond and silicon at various 
levels of theory (cf. text). 





diagonal 


off-diag 


;onal 




diamond silicon 


diamond 


silicon 


Hartree-Fock 


14.77 7.07 


0.23 


0.20 


relaxation / polarization 


-1.28 -0.98 


-0.01 


-0.02 


loss of ground-state corr. 


0.20 0.13 






long-range rel./pol." 


-3.5 -2.4 






a estimate based on the date in Ref. 25 
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(5.7 for diamond and 12.0 for silicon) 42 the above esti- 
mates of the long-range polarization contributions corre- 
spond to cutoff radii of 1.72 A for diamond and 2.75 A 
for silicon. To gain some insight into the accuracy of the 
long-range estimates these cutoff radii can be compared 
to the radii of the 8 atomic cluster kernels (defined here 
as the distance of the outermost kernel atoms from the 
center of the kernel). They amount to 1.89 and 2.88 A 
for C and Si, respectively, which is fairly close. 

Turning now to the off-diagonal matrix elements be- 
tween nearest-neighbor bonds we find again here that 
the correlation-induced changes are much smaller than 
for the corresponding valence band quantities evaluated 
in Refs. [ 14, 24] (changes in the matrix elements of 0.7 
cV for C and 0.8 eV for Si). In diamond lattices the 
nearest-neighbor off-diagonal elements enter with a factor 
of 8 into the bandwidth at the T-point of the uppermost 
(lowermost) four valence (conduction) bands and with a 
factor of 2 into the energy at the top (bottom) of the 
respective band complex. It is essentially these contribu- 
tion from the nearest-neighbor matrix elements which is 
responsible for the total correlation-induced change of the 
valence-band width (of ^7 cV for both, C and Si 14,24 ). 
Assuming a similar finding for the conduction band, i.e., 
that only matrix elements between nearest neighbors 
are important we would estimate the correlation-induced 
narrowing of its width to be in the order of 0.1-0.2 eV. 
The rising of its bottom would remain well below 0.05 
eV, an almost negligible shift. Yet, this assumption is 
somewhat unrealistic, because the larger extent of the 
conduction band WFs will lead to a much longer range 
of the correlation-induced changes of the local matrix el- 
ements than in the valence band case, i.e., to substantial 
changes far beyond nearest neighbours. This is corrob- 
orated by our ongoing study of the conduction bands of 
diamond by means of the local Hamiltonian approach 34 . 
Unfortunately, at present we do not have estimated data 
for these changes available. 

Indeed, our estimates for the center-of-gravity shifts 
(cf . above) , when taken together with the top-of- valence- 
band shifts of Rcf. 25 (i.e., 3.90 eV for C and 2.63 cV for 
Si) would lead to significantly too small band gaps (5.3 
eV vs. 7.3 eV from experiment 43 for C, and 2.5 eV vs. 3.4 
eV from experiment (direct gap) 43 for Si). They are cer- 
tainly more accurate than the (too large) HF band gaps 
of 13.8 eV for C, and 8.4 eV for Si 44 , but are still not ac- 
ceptable. More accurate data for the correlation-induced 
narrowing of the conduction band width are needed, and 
also the determination of the polarization contributions 
beyond nearest-neighbor bonds must be improved. Work 
in these directions is underway in our laboratory. 
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